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We study the production of charmonia and charmed hadrons for nucleus-nucleus colli- 
sions at SPS, RHIC, and LHC energies within the framework of the statistical hadroniza- 
tion model. Results from this model are compared to the observed centrality dependence 
of J/xIj production at SPS energy. We further provide predictions for the centrality de- 
pendence of the production of open and hidden charm mesons at RHIC and LHC. 

1. Introduction 

The original idea of statistical hadronization of charm quarks has sparked an 
intense activity in this field Initial interest focussed on the available SPS data on 

J/ip production. As we show below, these data can be well described, but only assuming 
an enhanced charm cross section compared to pQCD calculations. However, the largest 
differences between results from the statistical coalescence scenario (or a similar model 
0) and more conventional models are expected at collider energies. For example, in the 
Satz-Matsui approach 0, one would expect very strong suppression compared to direct 
production of J/ip mesons (up to a factor of 20 for central Au-Au collision at RHIC 
energy. In the present approach this suppression is overcome by statistical recombination 
of J/ip mesons from different cc pairs, so that much larger yields are expected. We 
therefore focus in this note on predictions for open and hidden charm mesons at RHIC 
and LHC energy, with emphasis on the centrality dependence of rapidity densities. 



2. Model description 

The model assumes that all charm quarks are produced in primary hard collisions and 
equilibrate in quark-gluon plasma (QGP); in particular, no J/ip is preformed in QGP 
(complete screening). All charmed hadrons (open and hidden) are formed at freeze-out 
(at SPS and beyond, freeze-out is at the phase boundary [||) according to the statistical 
laws. Taking into account the measured dependence of the ratio on centrality, 

the model is valid for not too small number of participants, Npart- The total number of 
thermally produced open charm hadrons, A*"*^ is related to the number of directly produced 
cc pairs, N^i"^ as (neglecting charmonia): N^}^'^ = \gcN^J^Ii{gcN^J^) / Iq^QcNI^), from which 
the charm enhancement factor is extracted. J„ are the modified Bessel functions. The 
yield of a given species X is then Nx = gcNxhigcNx)/ Io{9cNx) for open charm mesons 
and hyperons and Nx = Oc^x charmonia (see ref. for more details). The inputs 
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for the above procedure are: i) the total charged particles yields (or rapidity densities), 
which are taken from experiments at SPS and RHIC and extrapolated for LHC; 

and ii) N^^"^ , which is taken from next-to-leading order (NLO) perturbative QCD (pQCD) 
calculations for pp [|T^ (the yield from MRST HO parton distributions was used here) and 
scaled to AA via the nuclear overlap function. A constant temperature of 170 MeV and the 
baryonic chemical potential fib according to parametrization fif, = 1270/(1 + ^/sN^/4:.3) 
[0 have been used for the calculations. 

3. Results 

We first compare predictions of the model to 47r-integrated J/ip data at the SPS (NA50 
data replotted as outlined in [0]). The differences between this calculation and that of 



[g] are due to the inclusion in the present version of the complete set of charmed mesons 
(open and hidden) and baryons as well as due to updated values of the inputs (volume and 
open charm cross section) as outlined above. In Fig. |l] we present the comparison for two 
values of N^^"^: from NLO calculations |T2| and scaled up by a factor of 2.8. Note that the 
observed centrality dependence is well reproduced using the NLO cross sections for charm 
production scaled by the nuclear overlap function. To explain the overall magnitude of 
the data a A^'J^ increase by a factor of 2.8 compared to NLO calculations is needed. 
We mention in this context that the observed enhancement of the dilepton yield at 



intermediate masses [|T5| has been interpreted as possible indication for an anomalous 
increase in the charm cross section. If, in our calculations, we scale-up A^^^ from NLO 



according to the NA50 enhancement ("Data/Expected sources") as a function of Npart [|15 
the J/tp data are underpredicted and the Npart dependence is flatter (see Fig. p. However, 
the NA50 suggestion of charm enhancement is controversial, and other explanations exist 



of the observed enhancement in terms of thermal radiation 16,17 
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Figure 1. The centrality de- 
pendence of J/ip production 
at SPS. Model predictions are 
compared to NA50 data (47r- 
integrated). Two curves for 
the model correspond to values 
of N^^"^ from NLO calculations 
[0] (0.137, leading to c/c=0.78) 
and scaled up by a factor of 2.8. 
The dotted curve is obtained 
when considering a charm en- 
hancement as suggested by 
NA50 |T§. For central colli- 
sions (Apa,.f=350) we use in the 
calculations V=3070 fm^ and 
A'*!!=0.98. 
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We turn now to predictions for collider energies. For comparison we include in this 
study also results at SPS energy. The input parameters for these calculations for central 
collisions (A'parj=350) are presented in Table ^ Notice that from now on we focus on 
rapidity densities, which are the relevant observables at the colliders. The results are 
compiled in Table |^ for a selection of hadrons with open and hidden charm. All predicted 
yields increase strongly with energy, reflecting the increasing charm cross section and the 
concomitant importance of statistical recombination. Also ratios of open charm hadrons 
evolve with increasing energy, reflecting the corresponding decrease in charm chemical 
potential. 



Table 1 

Input parameters for model calcula- 
tions at SPS, RHIC and LHC. 



Table 2 

Results of model calculations at SPS, RHIC and 

LHC for Nr,art = ^50. 

200 55(Xr 



^SNN (CeV) 


17.3 200 5500 


T (MeV) 
fib (MeV) 


170 170 170 
253 27 1 


dN,h/dy 
VAy=i (fm^) 


430 730 2000 
861 1663 4564 


dN^^^'/dy 
9c 


0.064 1.92 16.8 
1.86 8.33 23.2 



y/SNN (GeV) 



dNT)+/dy 
dNj,-/dy 
diVoo /dy 
dN^Jdy 
dNj/^/dy 
dN^'/dy 



17.3 



0.010 
0.016 
0.022 
0.014 
2.55-10-^ 
0.95-10-5 



0.404 3.56 

0.420 3.53 

0.89 7.80 

0.153 1.16 

0.011 0.226 

3.97-10-^ 8.46-10-3 




Figure 2. Centrality dependence of J/ip ra- Figure 3. Centrality dependence of D'*' ra- 
pidity density at SPS, RHIC and LHC. pidity density at SPS, RHIC and LHC. 



Predictions for the centrality dependence of J/ ip production are presented in Fig. ^. In 
addition to the dramatic change in magnitude (note the scale-up by factors of 10 and 100 
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for RHIC and SPS energy, respectively) the results exhibit a striking change in centrality 
dependence, reflecting the transition from a canonical to a grand-canonical regime (see 
[0] for more details). The preliminary PHENIX results on J/ip production at RHIC []T8| 
agree, within the still large error bars, with our predictions. A stringent test of the 
present model can only be made when high statistics J/ip data are available. Another 
important issue in this respect is the accuracy of the charm cross section, for which so 
far only indirect measurements are available [0. In any case, very large suppression 
factors as predicted, e.g., by seem not supported by the data. In Fig. ^we present the 
predicted centrality dependence of charged D"'"-meson production for the three energies. 
The expected approximate scaling of the ratio D'^ /Npart like Np^^f. (dashed hues in Fig. ||) 
is only roughly fulfilled due to departures of the nuclear overlap function from the simple 
Nplrt dependence. 

4. Conclusions 

We have demonstrated that the statistical coalescence approach yields a good descrip- 
tion of the measured centrality dependence of J/ip production at SPS energy, albeit with 
a charm cross section increased by a factor of 2.8 compared to current NLO calculation. 
Rapidity densities for open and hidden charm mesons are predicted to increase strongly 
with energy, with striking changes in centrality dependence. First RHIC data on J/ip pro- 
duction support the current predictions, although the errors are too large to make firm 
conclusions. The statistical coalescence implies travel of charm quarks over significant 
distances in QGP. If the model predictions will describe consistently precision data this 
would be a clear signal for the presence of a deconfined phase. 
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